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Abstract. Preliminary results are presented about a fully self-consistent ,/V-body simulation 
of a sample of four massive globular clusters in close interaction within the central region of 
a galaxy. The Af-body representation (with N = 1.5 x 10 6 particles in total) of both the clus- 
ters and the galaxy allows to include in a natural and self-consistent way dynamical friction 
and tidal interactions. The results confirm the decay and merging of globulars as a viable 
scenario for the formation/accretion of compact nuclear clusters. Specifically: i) the fric- 
tional orbital decay is ~ 2 times faster than that predicted by the generalized Chandrasekhar 
formula; ii) the progenitor clusters merge in less than 20 galactic core-crossing times fe); 
iii) the NC configuration keeps quasi-stable at least within ~ 70t b . 
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1. Introduction 

Nuclear clusters ar e common acro s s the 
Hubble sequence dMatthews et al. I Il999t 
ISchodel et all l2007l) . In particular, th e ACS 
Virgo Cluster Survey dCote et al. I |2006|) 
shows many compact nuclei at the pho- 
tocenters of many of early-type galaxies 
(> 66% for My < -15). Their half-mass 
radii (r/,) are in the 2 4- 62 pc range, with 
average value (r;,) = 4.2 pc, and scale with 
the nucleus luminosity as r/, oc £0-5±o.oo3_ 
Brighter (My < -20.5) core-Sersic galaxies 
lack resolved stellar nuclei. The mean of the 
frequency function for the nucleus-to-galaxy 
luminosity ratio in nucleated galaxies, Log 
rj = -2.49 + 0.09 is indistinguishable from 
that of the Super massive black hole-to-bulge 
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mass ratio, Log (M./M ga i) = -2.61 + 0.07, 
calculated in 23 early-type galaxies with 
detected supermassive black holes (SBHs).It 
is thus argued that resolved stellar nuclei are 
the low-mass counterparts of nuclei hosting 
SBHs detected in the bright galaxies. If this 
view is correct, then one should think in terms 
of central massive objects, either SBHs or 
compact stellar clusters (CSCs), that accom- 
pany the formation and/or early evolution of 
almost all early-type galaxies. Comparing 
the nuclei to the nuclear clustersof late-type 
spiral galaxies reveals a close match in terms 
of size, luminosity, and overall frequency. A 
formation mechanism that is rather insensitive 
to the detailed properties of the host galaxy 
properties is required to explain this ubiquity 
and homogeneity. 

Another observational relevant point is the 
growing evidence of presence of very massive 
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(> 10 7 M ) YOUNG star clusters in Antennae, 
Magellanic Clouds, M 3 3, M8 2 Fornax dSph 
(Fritze-v. Alvensleben 1999; Ide Griis et al. I 
l2005t iFusi Pecci et al. Il2005h . as well as OLD 
dHarris & Pudritzlll994 in M87 and Virgo el- 
lipticals; 

lHarris et al, ] (120061) indicate how up to 
a 40% of the total mass in GCS of bright- 
est cluster galaxies is contributed by mas- 
sive (p. d. mass > 1.5 x 10 6 M G , in good 
agreement with re c ent th eoretical results by 
iKratsov & Gnedin I (120051) . Putting together 
these observational data (presence of resolved 
stellar nuclei in galaxies and likely initial pres- 
ence of massive star clusters, raise two ques- 
tions and possible answers. The questions are: 
i) how are these 'stellar' nuclei formed?, and ii) 
why very massive stellar clusters are no more 
observed in many galaxies? The (possible) an- 
swers could be resumed in: resolved stellar 
nuclei are formed via massive stellar clusters 
merging, after a substantial orbital decay to- 
ward the galactic central region. In the follow- 
ing we describe some of our work done re- 
cently to give substance to this explanation. 

2. Nuclei as remains of merged 
Globular Clusters? 

Many papers dealt with dynamical friction on 
massive objects orbiting galaxies. It seems well 
ascertained that sufficiently massive and com- 
pact clusters may decay towards their parent 
galaxy central region in a time short respect to 
the Hubble time. The less symmetries in the 
galactci potential, the quicker the orbital de- 
cay, tha depends (obviously) also on the initial 
orbital energy and angular momentum. An en- 
hancement of the classical frictional decelera- 
tion caused by background stars is due to tidal 
braking torque caused by the galactic potential. 
As a consequence, it has been proved in quite 
general cases that after less than 1 Gyr many 
GCs are limited to move in the inner galactic 
region. Do they merge and form a Compact 
Star Cluster? The first positive answer to this 
question was given by Trem aine et al. I d 19751) 
who were the first to examine the modes 
of cluster merging in the inner region of M 
31. Here, we summarize the results presented 



in ICapuzzo-Dolcetta & Miocchi (2008b) of a 
fully self-consistent iV-body simulation of the 
close interaction of a sample of four massive 
globular clusters (GCs) in the central region 
of a galaxy. Both the clusters and the galaxy 
are represented by mutually interacting par- 
ticles, thus including in a natural and self- 
consistent way dynamical friction and tidal 
interactions. This study represents a substan- 
tial improvement in the analysis of the fric- 
tional decaying and merging of GCs in galac- 
tic nuclear regions, a scena rio first tackled 
by semi-analitical approaches ( Tremaine et ai~l 
ll975UCapuzzo -Dolcettalll993h an d then pur- 
sued by jV-body experiments dOh & Lin 
l2000t ICapuzzo-Dolcetta & Miocchi I l2008al) . 
Clarifying the role of the above-mentioned dy- 
namical effects is important also to understand 
the formation and origin of Nu c lear Clusters 
(NCs ) (e.g. lOh&Linl 120001: iBekki et al. I 
12004 . 

3. Models and Results 

Each GC is represented by 256,000 particles 
initially distributed according to a King pro- 
file whose structural parameters are taken from 
th e set of the most compact clusters sim ulated 
in lCapuzzo-Dolcetta & Miocchi 1 d2008 al). The 
GCs are initially located at rest within 
the galactic core (see Fig. Q]). The galac- 
tic model is given by a spherical and 
isotropic Plummer phase-space distribution 
sampled with 512,000 particles. The simula- 
tion is performed with our own parallel tree- 
code using individual and variable time-steps 
dMiocchi & Capuzzo-Dolcettall2002l) . 

The simulation results can be re-scaled 
with any given set of galactic structural param- 
eters. One possible choice for these parameters 
is the following: core radius rt, = 200 pc; core- 
crossing time tb = 0.54 Myr; central density 
pto = 370M o pc- 3 . 

The main results of the simulation can be 
summarized as follows: i) the frictional orbital 
decay is ~ 2 times faster than that given by the 
use of the generalized Chandrasekhar formula; 
ii) the progenitor clusters (initially located 
within the galactic core) merge in less than 
20 galactic core-crossing time (~ 11 Myr), 
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Fig. 1. Evolution of the Lagrangian radii of the four Fig. 2. Upper panel: time evolution of the distance 
GCs as a whole. of the cluster CD from the galactic centre (rCD). 



see Fig . Q] and ICapuzzo-Dolcetta & Miocchil 
(2008b); iii) the NC configuration is quasi- 
stable at least within the simulated time 
(~ 70*6 ~ 40 Myr); iv) the total sur- 
face density profile has the typical appear- 
ance of a nucleated galaxy central profile, 
see Fig. [3j v) the global velocity disper- 
sion profile decreases towards the centre as 
found in the iGeha et ai~l d2002l) observations. 
These results are described in more detail in 
Cap uzzo-Dolc etta & Mi occhi 1 d2008bl) . 
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(solid line) and, for the sake of comparison, of the 
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GC progenitors. 
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Fig. 4. Surface brightness (g-band) profile of the 
dwarf ellipt ical galaxy VCC 1871 in the Virgo clus- 
ter (from JCote et al. l2"00g) 1. 
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Fig. 5. Projected radial behaviour of the velocity 
dispersion in the last system configuration. Line 
symbols are as in Fig. 4. 
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